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ABSTRACT 

Recent MHD dynamo simulations for magnetic Prandtl number > 1 
demonstrate that when MHD turbulence is forced with sufficient kinetic helicity, 
the saturated magnetic energy spectrum evolves from having a single peak 
below the forcing scale to become doubly peaked with one peak at the system 
(=largest) scale and one at the forcing scale. The system scale field growth is 
well modeled by a recent nonlinear two-scale nonlinear helical dynamo theory 
in which the system and forcing scales carry magnetic helicity of opposite sign. 
But a two-scale theory cannot model the shift of the small-scale peak toward 
the forcing scale. Here I develop a four-scale helical dynamo theory which shows 
that the small-scale helical magnetic energy first saturates at very small scales, 
but then successively saturates at larger values at larger scales, eventually 
becoming dominated by the forcing scale. The transfer of the small scale peak 
to the forcing scale is completed by the end of the kinematic growth regime 
of the large scale field, and does not depend on magnetic Reynolds number 
Rm for large Rm- The four-scale and two-scale theories subsequently evolve 
almost identically, and both show significant field growth on the system and 
forcing scales that is independent of Rm- In the present approach, the helical 
and nonhelical parts of the spectrum are largely decoupled. Implications for 
fractionally helical turbulence are discussed. 

Key Words: MHD; turbulence; ISM: magnetic fields; galaxies: magnetic fields; stars: 
magnetic fields; methods: numerical 



1. Introduction 



Most all ionized astrophysical rotators are turbulent and carry magnetic fields. For 
astrophysical objects such as the sun, helical dynamo theory was original proposed as a 
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theory to model the large scale field generation on scales larger than that of the driving 
turbulence (e.g. Parker 1955; Steenbeck, Krause, Radler 1966; Moffatt 1978; Parker 1979). 
Also of interest is to understand the overall shape and amplitude of the magnetic energy- 
spectrum from helical or nonhelical dynamo action on both large and small scales (e.g. 
Pouquet, Frisch, Leorat 1976, Kulsrud & Anderson 1992, Brandenburg et al. 1995; Hawley, 
Gammie, Balbus 1995,1996; Brandenburg 2001, Shekochihin et al. 2002ab, Maron & 
Cowley 2002; Maron & Blackman 2002). When the effect of the growing magnetic field on 
the turbulent velocity is ignored, the dynamo model is kinematic. We know that kinematic 
theory is incomplete since the dynamo must represent a solution of the full nonlinear 
magnetohydrodynamic (MHD) equations. In this context, much of theoretical dynamo 
research can be divided into two categories: (1) that which focuses on the fundamental 
principles, with the goal of determining a correct and useful nonlinear theory, and (2) 
that which parameterizes the backreaction of the growing magnetic field on the turbulent 
velocities and uses it in the dynamo equations to model the fields of specific objects in 
detail. This paper falls into the former category. 

Turbulence that produces astrophysical dynamos can either arise from external forcing 
or be self-driven from shear. Examples of the former are supernova driving in the Galaxy 
(e.g. Ruzmaikin, Shukurov, Sokoloff 1988) or thermal convection driving in the sun (e.g. 
Parker 1979), and an example of the latter is the magneto-rotational instability (MRI) 
(Balbus & Hawley 1991; 1998). In either forced or self-driven dynamos, there exist both 
helical or nonhelical varieties. Nonhelical vs. helical refers to the situation in which 
turbulence has either zero vs. finite pseudoscalar quantities such as the kinetic helicity 
v ■ Vxv, where v is the turbulent velocity, and the overbar represents some volume average 
(within a given hemisphere for a rotator). Such a correlation can be derived from the 
statistical tendency for rising (falling) eddies to twist in the opposite (same) sense as the 
underlying rotation as they expand (contract) and conserve angular momentum in stratified 
rotator. When kinetic helicity is present, the turbulence can generate a global large scale 
field (e.g. Moffatt 1978; Parker 1979) whose associated sign of flux lasts many overturn 
times of the largest turbulent eddy. As discussed in Blackman (2002b), such a global field 
is not to be confused with fields produced in nonhelical simulations for which the largest 
turbulent scale is the box scale. In this latter case, magnetic energy will arise on the scale 
of the box, but the sign of its flux will change every overturn time scale. For example, 
nonstratified MRI simulations (Hawley, Gammie, Balbus 1995) for which no kinetic helicity 
is present, produce magnetic energy with maximum azimuthal scale of order the box height, 
but the turbulence also extends to this scale, and the field does not maintain a coherent 
flux longer than a large eddy turnover time (~ an orbit time). In contrast, helical MRI 
simulations in a non-periodic box (Brandenburg et al. 1995) were shown to incur a helical 
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dynamo which sustains a large scale flux over many orbit times. This latter type of global 
field is extremely helpful in producing jets and coronae, since it can survive the buoyant rise 
without shredding on its way up (Blackman 2003). (Note that Stone et al. (1996) claimed 
that no helicity was observed in stratified periodic box disk simulations. But they averaged 
over their entire box to test for helicity, and since the northen and southern hemispheres 
must produce opposite sign, the box averaged helicity would be expected to vanish.) 

As the amplification of magnetic fields and the associated MHD turbulence represent 
highly nonlinear problems, understanding the backreaction of the magnetic field on the 
field growth itself has been an evolving topic of study (e.g. Piddington 1981, Kleeorin 
& Ruzmaikin 1982, Vainshtein & Cattaneo 1992, Cattaneo & Hughes 1996, Kulsrud 
& Anderson 1992, Field, Blackman, Chou 1999, Blackman & Field 2000, Brandenburg 

2001, Field & Blackman 2002, Blackman & Field 2002). Some recent developments 
(Brandenburg 2001; Field & Blackman 2002; Blackman & Field 2002; Brandenburg, 
Dobler, Subramanian 2002; Blackman & Brandenburg 2002; Maron & Blackman 2002) have 
emerged from focusing numerically and analytically on the simplest dynamo forced dynamos 
for which the underlying principles can be identified. Experiments have been performed in 
which a turbulent velocity spectrum is first established in the absence of magnetic fields in 
a periodic box and then a weak seed spectrum of magnetic energy is input (e.g. Meneguzzi, 
Frisch, Pouquet 1981; Kida, Yanase, Mizushima 1991; Brandenburg 2001, Maron & Cowley 

2002, Maron & Blackman 2002). The nonlinear evolution the magnetic energy spectrum to 
saturation can then be studied as the system is driven. (Note that these simulations are 
distinct from those in which an initial uniform background field is also assumed e.g. Maron 
& Goldreich 2002; Cho, Lazarian & Vishniac 2002) 

The shape of the resulting spectrum depends on the fractional helicity. Define the 
fractional kinetic helicity f h = |(v/ • Vxv/) /vjk 2 \ : where vj is the turbulent velocity on the 
forcing scale and kf is the forcing wavenumber. When forced such that fh = 0, the magnetic 
spectrum saturates with the magnetic energy seemingly piling up close to the resistive scale 
for magnetic Prandtl number (the ratio of viscosity to magnetic diffusivity Ptm = v/X) 
> 1, not near the forcing scale (Kida, Yanase, Mizushima 1991; Maron & Cowley 2001; 
Schekochihin 2002ab). This contradicts, for example, observations of the Galactic magnetic 
field, which does not seem to have a peak on the resistive scale, (Beck et al. & 1996). 
On the other hand, recent simulations of Haugen et al. (2003) suggest that the peak may 
not be at the resistive scale even for Pr > 1. Determing the location of the peak for the 
nonhelical dynamo for Pr > 1 is presently an active area of research. All studies seem to at 
least agree that the peak is below the forcing scale. 

By contrast, simulations that force with fh ~ 1 show production of large scale fields at 
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the box scale and a peak at the forcing scale (Brandenburg 2001). Maron and Blackman 
(2002) studied what happens to the spectrum both above and below the forcing scale as a 
function of fh- By starting with an initial spectrum that represented the saturated state 
of an fh = simulation, they found that for fh above a critical value (= kxjki, where 
k\ is the box scale) a peak at k\ emerged, and as fh was further increased toward 1, the 
peak originally at the sub-forcing scale increasingly depleted. A peak at the forcing scale 
emerged. The adjusted shape of the spectrum at and below the forcing scale established 
itself on a kinematic growth time scale of the k\ field. 

The growth rate and saturation value of the large scale field at k\ is now well 
understood by a nonlinear two-scale dynamo theory based on magnetic helicity conservation 
and exchange between the small scale (assumed to be the forcing scale) and the box (large) 
scale (Field & Blackman 2002; Blackman & Field 2002): Growth of the large scale field 
corresponds to a growth of large scale magnetic helicity. Magnetic helicity conservation 
dictates that a small scale helical field must then also grow of the opposite sign. Because 
the growth driver for the large scale helicity also depends on the small scale magnetic 
helicity, the growth of the latter ends up quenching the large scale dynamo. But there is 
an important limitation of the two-scale theory. Why should the growth of the small scale 
compensating magnetic helicity occur at the forcing scale and not at some much smaller 
scale? Answering this can help to understand why the small-scale magnetic energy peak 
initially at sub-forcing scales in simulations of Maron and Blackman (2002) migrates to the 
forcing scale when sufficient kinetic helicity is input. 

To address these questions, I develop a simple four-scale approach which predicts that 
the small scale magnetic helicity and current helicity grow fastest on the smallest scale 
on which there is both a finite turbulent velocity and magnetic field. (This would be the 
viscous scale for Pr M > 1 and the resistive scale for Pr M < 1.) The magnetic engergy and 
current helicities saturate fastest at these small scales. Eventually however, larger values 
are reached at larger scales until the small scale growth is dominated by the forcing scale, 
essentially justifying a two-scale approach by the time the kinematic growth phase for the 
k\ field ends. It is important to emphasize that the mean field formalism herein describes 
the dynamical evolution of both the large and small scale magnetic and current helicities. 
The backreaction of the small scale current helicity on the kinetic helicity, and thus on the 
growth of the large scale field is treated dynamically. We will see that significant growth of 
the large scale field proceeds unimpeded at early times, but at late times, the backreaction 
from the small scale field slows the growth to a resistively limited pace. It is imporant to 
emphasize that the dynamical mean field theory studied herein, albeit a simplified theory, 
goes beyond the kinematic mean field theory of standard textbooks which does not include 
the backreaction. 
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In section 2, the basic role of helicity conservation in large scale dynamo theory is 
discussed and the four-scale set of equations (3 of which are coupled) are derived. In section 
3, the equations are solved. The solutions and their interpretation, discussed in section 4. 
Section 5 is the conclusion. 



2. The Four-scale Nonlinear Dynamo 

The mean field dynamo (MFD) theory has been a useful framework for modeling the in 
situ origin of large-scale magnetic field growth in planets, stars, and galaxies (Moffatt 1978; 
Parker 1979 Krause & Radler 1980; Zeldovich, Ruzmaikin, Sokoloff 1983) and has also been 
invoked and to explain the sustenance of fields in fusion devices (see Ortolani & Schnack 
1993; Bellan 2000; Ji & Prager 2002 for reviews). Pouquet, Frisch & Leorat (1976) derived 
approximate evolution equations for the spectra of kinetic energy, magnetic energy, kinetic 
helicity, and magnetic helicity and demonstrated the need to consider magnetic helicity 
in more carefully in the context of dynamo theory. Recent progress in modeling the fully 
nonlinear dynamical field growth and saturation seen in simulations (Brandenburg 2001) 
has come from dynamically incorporating magnetic helicity evolution using a two-scale 
theory (Field & Blackman 2002, Blackman & Brandenburg 2002, Blackman & Field 
2002). Previous attempts (Zeldovich, Ruzmaikin, Sokoloff 1983, Kleeorin & Ruzmaikin 
1982, Kleeorin, Rogachevskii, Ruzmaikin 1995, Gruzinov & Diamond 1995, Bhattacharjee 
& Yuan 1995) to use magnetic helicity in dynamo theory did not produce a time dependent 
dynamical theory. 

To develop the new four scale theory and make contact with the previous work, I first 
average the magnetic induction equation to obtain the basic MFD equation (Moffatt 1978) 

d t B = Vx£ + Vx(V x B) - AV 2 B, (1) 

2 

where B is the mean (large-scale) magnetic field in Alfven speed units, A = ^ is the 
magnetic diffusivity in terms of the resistivity 77, V is the mean velocity which I set = 0, 
and £ = (v x b) is the turbulent electromotive force, a correlation between fluctuating 
velocity v and magnetic field b in Alfven units. Instead of the assuming that there are 
two scales as in Field & Blackman (2002) and Blackman & Field (2002), here I generalize 
and assume that there are four scales on which magnetic helicity can reside, such that 
&4 > k 3 > k 2 > ki > 0. An example is to take k 2 = kf, the forcing wavenumber, fc 4 = k\, 
the resistive wavenumber, and k^ = k u viscous wavenumber. We can interpret B, A and £ 
as the mean components of B, A and (v x b) that varies on a scale k^ 1 in a closed system. 

The total magnetic helicity, H = (A-~B) vol , satisfies (Woltjer 1958; Moffatt 1978 Berger 
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& Field 1984) 

d t H = -2(E • B) voh (2) 

where E = — <9 t A — V0, is the scalar potential, and () vo i indicates a global volume 
average. Such a conservation formula also applies separately to each scale. That is, using 
E = <9fA — V0, e 2 = d t SL 2 — V0 2 and e 3 = d t a 3 — V0 3 and dotting with B, b 2 and b 3 
respectively we obtain 

d t E x = 2(£- B) vol - 2A(J • B) vo i, (3) 
d t H 2 = -2(e 2 • b 2 )„ oi , (4) 
d t H 3 = -2(e 3 • b 3 )„ oi , (5) 

and 

d t H 4 = -2(e 4 • b 4 )„ oi , (6) 

where ffi = (A • B)„ oi , H 2 = (a 2 ■ b 2 ) t)oi , # 3 = (a 3 • b 3 ) vol and H 4 = (a 4 • b 4 )„ oi . I now 
derive expressions for {£ ■ B)„ oJ , (e 2 • b 2 )^, (e 3 ■ b 3 )„ j, and (e 4 • b 4 )^. I assume that mixed 
correlations between fluctuating quantities of widely separated scales (which fluctuate on 
widely different time scales) do not contribute. Thus I ignore terms of the form e.g. (b 2 • j 3 ). 
This is an important simplifying assumption which really needs to be tested numerically, 
and I will come back to it again in section 4. Note however that this assumption is not 
the same as ignoring the effect of the small scale fields on the large scale field growth, or 
ignoring the nonlinear backreaction altogether: Eqn. (3-6) along with (8-13) show that 
nonlinear couplings between disparate scales also occur because the time evolution for 
the large scale magnetic helicity depends non-linearly on the time evolution of the small 
scale helicities. These latter couplings are included, even when mixed correlations are not 
included. We then have 



(£ ■ B) vol = (S 2 ■ B) vol + (£ 3 ■ B) vol (7) 
where £ 2 = v 2 x b 2 and £ 3 = v 3 x b 3 , and 

-(e 2 • b 2 ) vol = -(£ 2 ■ B) vol + ((v 3 x b 3 ) 2 • b 2 ) vol - A(j 2 • b 2 ) vol (8) 

-(e 3 • b 3 ) vot = -((v 3 x b 3 ) 2 • b 2 )„ i - (£ 3 ■ B) vol - A(j 3 • b 3 ) vol (9) 

and 

(e 4 • b 4 )„ oi = A(j 4 • b 4 )„ oi (10) 

where the overbar indicates averages which vary on scales /cf 1 wavenumber, the () 2 indicates 
averages which vary on scales > h 2 x . Eqn. (10) follows because v 4 = 0, as I have assumed 
k4 is taken at the resistive scale and that mixed scale correlations vanish. 
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In what follows I will solve for the nonlinear time evolution of the magnetic helicity at 
each of the four scales. Blackman & Field (2002) showed that a proper nonlinear theory 
must technically include the time evolution of £ 2 and by generalization here, also the time 
evolution for £ 3 into the theory. However, they also showed that for the specific case 
for which triple correlations in d t S 2 are treated as damping terms with a damping time 
t 2 = l/k 2 v 2 , neglecting the time evolution of the turbulent electromotive force does not 
qualitatively influence the solution. Adopting this specific case, the implication of their 
result for the present theory is that 



£2 = j(h b 2 - v 2 • Vxv 2 )B - (3 2 3, (11) 

where (5 2 is a diffusivity computed from velocities at the k 2 scale. (The form of the 
diffusivity is not essential for the present discussion, and I will later scale it to its kinematic 
value.) Similarly, we also have 

r 3, 



£ 3 = v 3 xb 3 = ^(j 3 -b 3 )B-/3 3 J, (12) 
and 

(v 3 x b 3 ) 2 = — (j 3 • b 3 ) 2 b 2 - (/3 3 ) 2 j 2 , (13) 

where r 3 = l/k 3 v 3 , and I assume that there are no kinetic helicity contributions at the k 3 
scale. This is reasonable because I am focusing on the case in which the kinetic helicity is 
peaked at the forcing scale, and kinetic helicity does not cascade efficiently. If the k\ field 
is maximally helical, the current helicity, magnetic helicity and magnetic energy for the k\ 
scale are simply related by powers of k±. (However, the third term of (8) and the second 
term of (9) involve a subtly that will be explained below.) Combining equations (6-13) we 
have 

dtH, = 2^(k 2 2 H 2 + f h k 2 vl)k x H x + 2^(A; 1 # 1 A; 3 2 iJ 3 ) - 2(A + 2 + J 3 )k\H x , (14) 

d t H 2 = -2^{k 2 2 H 2 + UhvDhH, - 2^(k 2 H 2 k 2 3 H 3 f u ) + i$ 2 k\H x - 2g u ^ 3 ) 2 k 2 2 H 2 - 2Xk 2 2 H 2 , 

(15) 

d t H 3 = 2^(k 2 H 2 k 2 3 H 3 f u ) - 2^(k 1 H 1 k 2 3 H 3 ) + 2{3 3 k\H x + 2g u {{3 3 ) 2 k 2 2 H 2 - 2\k 2 3 H 3 , (16) 
and 

d t H A = -2\k\H A . (17) 

The last equation is decoupled from the others, and represents decay. Thus in present 
approximation scheme, where correlations between different scales vanish and k A = k\> k 3 
there is no helicity exchange with the resistive scale. Equation (17) can then be subsequently 
ignored. 
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The third and fifth terms of (15) and the second and fifth terms of (16) involve the 
quantities f u and g u . These are positive quantities and they come from the subtlety in dotting 
the second and third terms of (13) with b 2 : dotting the third term of (13) with b 2 gives the 
term oc ((/3 3 ) 2 j 2 • h 2 ) vo i = (fcjvoiij^^h ■b 2 ) vo i = (Ps)voi = (j§^h -b 2 ) vo i = (As) voi9uk 2 H 2 , 
where g u accounts for the deviation of jj^ 2 ^ from unity and I also assume that (3 2 = {(3 2 ) vo i 
and f3 3 = (flz) vo i- Similarly dotting the second term of (13) with b 2 and averaging gives a 
term oc ((j 3 • b 3 ) 2 hl) vol = (j 3 • h 3 ) vol (^^-b 2 2 ) vol = -(j 3 • b 3 ) vol f u k 2 H 2 so that f u accounts 

for both the fractional helicity at the k 2 scale and the deviation of the ratio fa, b ^ 2 from 
unity. Note that unlike b\, b\ does not enter these equations so that even if there is a large 
pile-up in magnetic energy at k 3 initially only its force-free component couples into the 
above equations. 



To write equations (14), (15), and (16) in dimensionless form, I define hi = Hi(k 2 /v. 
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h 2 = H 2 (k 2 /vl), h 3 = H 3 (k 2 /vj), R M = v 2 /\k 2 , r = tv 2 k 2 , use v 2 • Vxv 2 = -f h k 2 v%. I 
assume that f3 2 = q(hi)v 2 /3k 2 and j3 3 = q(hi)v 3 /3k 3 , where q{h\) is a quenching function 
for the diffusivity that I take to be q(hi) = (1 + kihi/k^ 1 . The main points of the present 
study are insensitive to the form of q, but pinning it down requires more study. I also use a 
Kolmogorov velocity spectrum so that v 3 = (k 2 /k 3 ) 1 ^ 3 v 2 and r 3 = {k 2 /k 3 ) 2 ^T 2 . The result is 

drh = l(h + h 2+ (|) 1 h 3 ) % hl - 2 + ^1 (I) 1 + ^) (|) 2 hl , (18) 

9rh = f (%ffuhsh 2 -§(/, + h 2 ) + » i (&) 2 ^l-2 (I) 1 + ^) h 2 , 

(19) 

and 



drh = I {&fuhh 2 - § (fc)* (|) h 3hl • ^-(fer (&) Z + (If) J h 2 

(20) 

These equations are to be solved for hi, h 2 and h 3 as a function of time. They can model 
a four-scale nonlinear dynamo theory in which helicity decays at the resistive scale but is 
allowed to transfer between the system scale, the forcing scale, and the viscous scale. We 
will compare the solutions with those of the two-scale system of equations: 

d T h 2 = --(f h + h 2 )-±h 1 + ^-^ -M h-—h 2 (22) 
3 k 2 o \k 2 J R M 
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obtained from (14), (15), and (16) by setting H 3 = and (3 3 = and then non- 
dimensionalizing. Equations (21) and (22) are essentially the same as those solved in 
Field & Blackman (2002) and which emerged as the specific limit of the strong damping 
approximation in the generalized two-scale theory of Blackman & Field (2002). Using the 
doublet (21) and (22) of equations rather than the triplet of equations which also includes 
d t £ cannot capture possible oscillations in the early time growth of hi. (Related oscillations 
were identifed in simulations by Stribling, Matthaeus, Ghosh 1994 when a uniform field 
was imposed in their periodic box.) These can be accounted for only when one allows for a 
variation in the damping time in the d£ 2 equation. Similarly, the use of the four equations 
(14-17) rather than 6 equations that would result from also including equations for d t £ 2 and 
d t £ 3 is similarly restrictive in the present context. We leave the necessary generalizations 
for future work. 



3. Discussion of Solutions 

The critical helicity for growth of hi to beat diffusion is fh > ki/k 2 . For f h below this 
value, hi, h 2 , and h 3 all decay. This is consistent with numerical simulations of Maron & 
Blackman (2002). Below we focus on the case fh — 1 to emphasize the main ideas. We also 
assume g u — f u — 1 to simplify the discussion, and a seed of hi ~ 0.001. The associated 
solutions of (18), (19), and (20) will be discussed along with the comparison to solutions of 
(21) and (22). 



3.1. Large scale field growth in the kinematic regime 

By kinematic regime we mean the growth regime in which hi grows independently of 
Rm- This is essentially the time for h 2 to approach —fh (see (18)). At early times, before 
resistivity is important, h 2 ~ —hi, so we can estimate the duration of the kinematic regime 
by determining the time at which hi ~ fh- We have, from (18), 

t kin ~3(k 2 /ki)Ln(f h /hi(0)), (23) 

which is ~ 100 for f h = 1, k 2 /ki = 5 and hi(0) = 0.001, as used herein. 

A careful look at (18), (19), and (20) shows that for a small seed hi and 
^3(0) %j fh(k 2 /k 3 )^ 3 , the kinetic helicity drives the growth of hi and the growth of the 
oppositely signed h 2 . (If h 3 (0) exceeds the above critical value then hi and h 3 initially 
decay until h 3 falls below that value, after which hi starts to grow.) The growth of h 2 also 
supplies, from diffusion, a growth of h 3 from the 5th term in (20). For k 2 < k 3 « k\, 
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the value of hi at the end of the kinematic phase in the four-scale approach saturates at a 
slightly smaller value of hi than for the two-scale case. The slight difference between the 
two curves can be explained by ignoring diffusion and dissipation terms (those containing q 
and Rm) in (18), (19), and (20) and assuming a quasi-steady state. Setting time derivatives 
to zero in (18), (19), and (20) then gives a system of three equations: 

(1 - h/k 2 + h 2 + (k 3 /k 2 ) 4 / 3 h 3 ) ~ 0, (24) 

h 2 h 3 = hih 3 (h/k 2 ) - h 2 (k 2 /k 3 f 3 ~ -h 2 {k 2 /k 3 f' 3 (25) 

and 

h 1 + h 2 + h 3 ~ 0, (26) 

where we have assumed h 2 » hiki/k 2 anticipating that h 3 is small. The solution 
is hi ~ 1 - h/k 2 - {k 2 /k 3 ) A ' 3 + (k 2 /k 3 f/ 3 , h 2 ~ -1 + (fc 2 /fc 3 ) 4/3 + h/k 2 , and 
/i 3 ~ —(k 2 /k 3 ) 8 / 3 . Thus no matter what k 3 « k\ is chosen, hi is only slightly 
depleted from hi ~ 1 — ki/k 2 , the value /ii attains in the kinematic regime for the two-scale 
approach. For k 3 — > A; 2 in (24), (25), and (26) /ii — > 1 — fci/fo. Also, for k 3 » 1, 
^! — > 1 — ki/k 2 . The correction to ft^ is maximized for ^2/^3 = 0.6 at which the deficit 
is about 25%. In short, hi — > 1, as in the kinematic regime here just as in the two-scale 
approach. For k 3 = 30 and Rm = 9000, these results are seen in the bottom pair of curves 
in the upper left plot of Fig. 1, where the bold line is the hi solution of the two-scale system 
(21) and (22) and the thin line is hi for the four-scale system. The top set of curves in the 
same plot shows the k 3 = 30 and Rm = 900 case. There, h 3 first grows as in the previous 
case, but drains more quickly from resistivity. Growth of hi then more closely mimics the 
kinematic evolution of hi in the two-scale approach at very early times. 

That h 3 ~ —(k 2 /k 3 ) 8 / 3 at the end of the kinematic regime shows that as k 3 increases 
for fixed k 2 , h 3 decreases. This is demonstrated in the bottom left plot of Fig. 1 which 
shows the growth of hi for R M = 100 and R M = 2 x 10 4 and k 3 = 160. There, even the 
high Rm four-scale solution matches the two-scale approach exactly. 

3.2. Rapid migration of the small scale peak to the forcing scale in the 

kinematic regime 

An important implication of including the scale k 3 can seen by comparing the ratio of 
the magnetic energy and current helicity at scales k 2 and k 3 to those of ki, as a function of 
time. This is shown in Figs. 2, 3 and 4 for different cases. In general, the current helicity 
and magnetic energy are dominated by the k 3 scale at early times. These quantities are 
related to the magnetic helicity by factors of k 2 and k respectively. The diffusion of h 2 into 
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h 3 grows enough h 3 for k 3 h 3 and k 3 h 3 to exceed the respective values at ki and k 2 initially. 
The early time dominance of the k 3 quantities is enhanced when there is an initial seed 
of h 3 , as seen by comparing Figs. 3 and 4. For a Kolmogorov velocity spectrum however, 
the small scale current helicity and associated magnetic energy end up being dominated by 
the forcing scale before tkm- as discussed below (26), the maximum h 3 at the end of the 
kinematic regime is h 3 ~ — (k 2 /k 3 ) 8 ^ 3 so that the ratio of current helicity at k 3 to that at k 2 
is then is {k 3 /k 2 ) 2 {h 3 /h 2 ) ~ {k 2 /k 3 ) 2 / 3 « 1, since h 2 ~ 1. 

It is interesting to asses whether the times for crossover of small scale dominance from 
k 3 to k 2 quantities depends on magnetic Reynolds number. The answer is that for large Rm 
it does not, but for low R M it does. One can be misled if applying the low Rm results to 
large Rm- To see this more explicitly, first note that the location of the crossovers for the 
case of Fig. 3 does depend on R M - The lower the Rm the earlier the crossovers. However, 
for large Rm, the location of the crossover asymptotes to an Rm independent value. It has 
been checked that the location of the crossovers for Rm = 9000 in Fig. 3 is indistinguishable 
for the location of the cross over for any higher Rm, keeping all other parameters fixed. 
This latter result is more directly seen in the comparisons of the top and bottom rows of 
Figs. 3 and 4. There, the top rows in both cases correspond to a comparison between 
Rm = 100 and Rm = 2 x 10 4 for k 3 = 160. The bottom rows correspond to a comparison 
between R M = 100 and Rm = 2 x 10 5 for k 3 = 160. The locations of the cross over for the 
R M = 2 x 10 4 and fi M = 2x 10 5 cases are indistinguishable. 

Note that the case of k 3 = 160 and Rm = 100 for k 2 = 5 is the case of k 3 ~ k\ ~ k u 
for Kolmogorov turbulence, in other words, Ptm — 1- This follows because choosing 
R M = 100 for k 3 = 160 also corresponds to choosing (k 3 /k 2 ) = i? 3 / 4 ~ 100, where R is 
the hydrodynamic Reynolds number and where k 3 is taken to be the viscous wavenumber. 
The power 3/4 follows from the Kolmogorov spectrum. This can be distinguished from any 
Rm » 100 for fixed k 3 , which, in turn, represents Ptm » 1- 

Thus, restating the results of this section in terms of Pr M , we can say that for the 
Pr M = 1 case the crossover time for the current helicity at k 3 to deplete below that at 
k 2 is earlier than for the Pr M » 1 case in which k 3 « k\. This is expected, since the 
resistivity is not effective at dissipating k 3 for Ptm >> 1- At large Rm for fixed k 3 = k v (or 
Ptm >> 1); the crossover is the result of the dynamical depletion of h 3 from the first two 
terms on the right of (20), and becomes independent of Rm or Pr M - 
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3.3. Saturation: the doubly maximal inverse transfer state 

The growth term for at early times is the fifth term in (20), but eventually this 
is offset by the fourth term, after which the only remaining terms for h 3 are resistive and 
inverse transfer loss terms. Thus, as we have seen, /i 3 eventually decays, whilst h 2 continues 
to grow. The latter takes over the role of compensating negative helicity, and the saturation 
proceeds exactly in the two-scale approach: positive helicity at hi and negative helicity h 2 
at the forcing scale. This is demonstrated in the right column of Fig. 1 where the late time 
evolution of the two-scale and four-scale approaches are seen to be indistinguishable for all 
Rm- Accordingly, at late times, the current helicity at k 2 asymptotes to equal that at k±, 
just as in the two-scale approach (Brandenburg 2001; Field & Blackman 2002; Blackman & 
Field 2002). 

The features of the four-scale model discussed herein can be used as an aid to 
understand the full saturated spectrum of the helical dynamo both above and below the 
forcing scale. Consider < k\ < kf < k 2 < k 3 < k u , where ki is the wavenumber of 
the largest scale available, kf is the forcing scale and k v is the viscous scale and we allow 
ki,k 2 , k 3 to take on intermediate values. The time for the spectrum to evolve depends on 
the ratios kj^/kf and kf/k v , but the qualitative evolution can be understood simply. At 
early times, a seed spectrum forced with negative kinetic helicity with f h > ki/kf, can 
grow positive magnetic helicity at k\. (For fh = 1 the maximum growth of the large scale 
magnetic helicity initially occurs at ki ~ k 2 /2.) The magnetic helicity will inverse transfer 
from k\ > k^ to successively smaller k, eventually all the way to k^. The ki will continue 
to accumulate the bulk of the positive magnetic helicity until saturation. This happens 
by direct analogy to the sub-forcing scale dynamics discussed earlier. There the negative 
magnetic helicity initially grows fastest on the scale k v before loss terms eventually win on 
that scale, and the negative magnetic helicity transfers upward to larger scales, eventually 
reaching the forcing scale where the bulk of the negative magnetic helicity collects. The 
positive and negative magnetic helicities thus migrate to their largest available scales. The 
negative helicity cannot grow at any k significantly less than kf since the growth rate 
for positive magnetic helicity is positive for lower k. The overall picture just outlined for 
saturation seems to be consistent with simulations of Brandenburg (2001) and Maron & 
Blackman (2002). 

Note that the larger the scale, the less important the resistive terms in draining 
the helicities from those scales. The dominant loss terms for k « k\ are the inverse 
transfer terms, not resistive terms. The threshold magnetic helicity at a given scale for the 
inverse transfer loss to dominate its gain from diffusion of magnetic helicity from larger 
scales, increases with decreasing k. Thus, larger scales take longer to inverse transfer their 
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magnetic helicity. For the positive helicity above the forcing scale, this means that the k L 
scale takes the longest to grow to its maximum, and for the negative helicity, kj takes the 
longest to grow to its maximum. The final state can be described as a state of "doubly 
maximal inverse transfer." 



3.4. Implications and comparison to fractionally helical dynamo spectra 

Maron & Blackman (2002) showed that as fh exceeds ki/kf and approaches 1 for 
Ptm > 1; the saturated magnetic energy spectrum changes from peaking on the resistive 
scale to peaking on two scales, with one peak at the forcing scale and one at the system 
scale. (Their fh = 1 results matched Brandenburg 2001.) For fh > ki/kf, fhiki/kf) times 
the equipartition energy ends up in magnetic energy at the the k\ scale after the kinematic 
regime. At saturation, fhk//k L of the equipartition energy ends up at k L . 

The fraction of small scale magnetic energy associated with the helicity dynamics is 
also determined by fh. Whatever the dynamics of the nonhelical magnetic energy spectrum, 
the helical fraction of the magnetic spectrum seems to be explicable by the independent 
dynamics discussed above, suggesting that the helical and nonhelical parts of the spectrum 
are rather decoupled. Recall that the nonhelical magnetic energies on the k 3 and k\ scales 
do not enter the theory, only that on the k 2 scale (see the discussion between equations (17) 
and (18)). In the present theory, the expected helical fraction of the small scale magnetic 
energy that follows the helicity dynamo dynamics is given simply by fh- This follows 
because the small scale current helicity in saturation is fhV 2 k 2 ~ (b • Vxb) and since 
(b 2 ) ~ (v 2 ) in saturation, we have the associated helical magnetic energy fraction = f\. If 
we include the nonhelical component of the magnetic energy, then fh is a lower limit to the 
fraction of the total small scale magnetic energy that winds up at the forcing scale. 

4. Conclusions 

A four scale nonlinear magnetic dynamo theory was presented using the approximation 
that correlations of mixed scales are assumed to vanish. The goal was to develop a simple 
theory that sheds light on the evolution of the full magnetic spectrum for the helical 
dynamo and show how the presence of kinetic helicity influences the spectrum both above 
and below the dominant scale of the turbulent kinetic energy. The velocity spectrum was 
assumed to be Kolmogorov, with kinetic helicity input only at the forcing scale. The results 
of the simple theory are consistent with existing numerical simulations and make additional 
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predictions which can be tested. The theory includes the dynamical backreaction of the 
growing magnetic field on the turbulence driving the field growth. 

The growth of the large scale field in the four scale theory is consistent with that 
predicted in the two-scale theory (Field & Blackman (2002); Brandenburg & Blackman 
(2002); Blackman & Field (2002)) at late times but previous simulations do not have 
enough resolution to test the dynamics of the kinematic regime and so higher resolution 
simulations will be needed. The four-scale theory herein predicts that as the large helical 
scale field grows, the small scale helical field of the opposite sign and its associated current 
and magnetic helicity will first grow at the smallest scales where both v and b are finite 
(for Pr>l, this is the viscous scale). This cannot be seen in a two-scale approach, since 
there the system scale and forcing scales are the only ones present. The source of helical 
field for the very smallest scales is diffusion from above. These very small (viscous) scale 
helical fields drain by inverse transfer, and the helical magnetic energy below the forcing 
scale peaks at successively larger scales, arriving at the forcing scale before tkm- If the 
viscous scale is much larger than the resistive scale, then resistivity plays little role in this 
process; the process is independent of Rm for large Rm- At tki n the growth of both the 
large and small scale helical fields proceed exactly as in the two-scale theory: large scale 
current and magnetic helicities are primarily carried by k±, and the oppositely signed small 
scale quantities at scale k 2 . 

That high resolution simulations are needed to test the kinematic regime also means 
that they are need for testing the assumption made above Eqn. (7). If future simulations 
for Ptm » 1 show that the current helicity piles up at the resistive scale rather than 
the viscous scale at very early times, this would suggest that the approximation that 
correlations of mixed scales vanish would have to be revised, since it is only such mixed 
correlations which can grow helicity on the resistive scale. In the present, four-scale theory, 
the helical magnetic energy on the resistive scale decays, and the nonhelical component is 
decoupled from the helical component. I emphasize that the model presented herein is a 
simplified theory, meant to guide physical understanding of helical MHD turbulence. It 
must be subject to further testing. 

The larger the fractional kinetic helicity fh, the more the overall magnetic spectrum 
behaves like its helical component. If the nonhelical component prefers to pile- up at small 
scales, the helical component will provide at least some fraction that wants to be doubly 
peaked at the forcing and system scales. It is thus important to emphasize that that driving 
turbulence with kinetic helicity on a scale k 2 affects the shape of the overall helical magnetic 
energy spectrum and thus that of the overall magnetic energy spectrum above and below 
the forcing scale. 
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Figure 1: The magnetic helicity, hi, for the four-scale approach (thin lines) and the 
two-scale approach (thick lines). The left column of plots is for early times and the right 
column is for a broader time range. For the top row, ki — 1, k 2 — 5, k 3 — 30, and Rm = 900 
(top pair of curves) Rm = 9000 (bottom pair of curves). For the bottom row, ki — 1, 
k 2 = 5, k 3 = 160, and Rm = 100 (top pair of curves) Rm = 2 x 10 4 (bottom pair of 
curves). Note that the two-scale and four-scale approaches are largely indistinguishable in 
all but the upper left plot (see text). (The high R M curve in the lower right will eventually 
saturate at the same value hi = (k 2 /ki) 2 as the low Rm case but at much later times and 
thus undershoots the top curve for the plotted time range.) 




Figure 2: The absolute magnitude of the current helicity ratios \k 2 h 2 /k\hi\ and \klh 2 /kfhi\, 
and helical magnetic energy ratios \k 2 h 2 /kihi\ and \k 2 h 2 /kihi\ assuming an initial seed of 
h 2 (0) = h 3 {0) = -h l {0)/2 = 0.0005. Here R M = 900 (thin lined curves) R M = 9000 (thick 
lined curves) for k\ — 1, k 2 — 5, k 3 = 30. In each plot, and for both values of Rm, the 
quantities at k 3 dominate at early times, and then become subdominant to the values at k 2 
at later times. The location of the crossover is earlier for Rm = 900 because the resistive 
scale is closer to the scale k 3 ; resistivity is more effective at early times in draining the 
current helicity and magnetic energy at k 3 than in the Rm = 9000 case. 
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Figure 3: The absolute magnitude of the current helicity ratios \k\hilk'\h\\ and | k\h,2/k\ h\\, 
and helical magnetic energy ratios ^2^2/^1 ^i| and | ^2/12/^1^1 1 assuming an initial seed of 
h 2 (0) = h 3 (0) = -/ii(0)/2 = 0.0005. for k x = 1, jfe 2 = 5, A; 3 = 160. Top row has plots for 
i?M — 100 (thin lined curves) and Rm = 2 x 10 4 (thick lined curves), and bottom row has 
plots for R M = 100 (thin lined curves) and Rm = 2 x 10 5 (thick lined curves). In each plot, 
the quantities at k% dominate at early times, but then become subdominant to the values at 
k 2 at later times. The magnetic energy plot is shown for a broader time range. The location 
of the crossover for the large R M cases of 2 x 10 4 and 2 x 10 5 is independent of R M - the 
crossovers occur for the thick pairs of lines at the same time in the top and bottom rows. 
The crossover occurs much earlier for Rm = 100 because the resistive wavenumber is closser 
to k% and thus resistivity is more effective at early times in draining the current helicity 
and magnetic energy at k% than in the larger Rm cases. For the current helicity, the cross 
over for Rm = 100 occurs so early that it is not visible on the graph. Before the end of the 
kinematic regime (t < 100), the crossovers are complete and fc 2 emerges as the dominant 
small scale. Since ki = 160 here, Rm = 100 corresponds to Pr M — 1 
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Figure 4: Same as Figure 3, but with ^(0) = 0.001 and h 2 (0) = h 3 (0) = 0. Note the initial 
rise of the current helicity and magnetic energy at ks but then again the dominance of the 
ki quantities before the end of the kinematic regime t^n = 100 



